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ABSTRACT 



We have reanalysed the giant outburst of the blazar PKS 15 10-089 (z 
the 0.1 - 100 GeV energy range exceeded the value of 10"' ph cm"^ ; 



: 0.36) occurred on 201 1 October-November The y-ray flux in 
' for several days. The peak flux was reached on 201 1 October 



variability was measured. Particularly, on 201 1 October 18, the flux doubling time was as small as ~ 20 minutes. This is the shortest 
variability ever detected in the MeV-GeV energy band. We compared the present analysis with two other outbursts observed in 2009 
March and 2012 February-March, when the blazar was detected also by HESS and MAGIC, respectively, in order to infer information 
about the emission at hundreds of GeV. 
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1. Introduction 

According to the current knowledge, the structure of relativis- 
tic jets in active galactic nuclei is self-similar (e.g. Heinz & 
Sunyaev 2003). This means that the size of the emitting region r 
is linked with the distance R from the central supermassive black 
hole with the relationship R ~ r/tfr, where iff is a scaling fac- 
tor generally with values between 0.1 and 0.25 (e.g. Ghisellini 
& Tavecchio 2009, Dermer et al. 2009), although radio obser- 
vations suggest even smaller values (e.g. Jorstad et al. 2005, 
Pushkarev et al. 2009). Since it is expected that the central black 
hole is the source of perturbation along the jet, the minimum 
size of the blob has to be larger than the gravitational radius 
rg = GM/c^, where M is the mass of the compact object, G 
the gravitational constant, and c is the speed of light in vac- 
uum. Given these constraints, the changes in the electromagnetic 
emission have to be on time scales t/(1 + z) > r/c6, where z is 
the redshift, and 3 is the Doppler factor of the jet. In the case of 
a typical blazar, with black hole mass of 10'' M©, Doppler factor 
6 ~ 10, and dissipation region located at R ~ lO^rg (Ghisellini 
et al. 2010) - which in turn means r ~ 2x lO'^ cm with tfr = 0.1 
- the observed time scale t has to be greater than ~ 18 hours. 

Shorter time scales have already been observed in the past 
with different facilities (Gaidos et al. 1996, Foschini et al. 2006, 
Aharonian et al. 2007, Albert et al. 2007). The launch of the 
Fermi Gamma-ray Space Telescope (hereafter Fermi) in 2008 
made it possible to increase significantly the sample of observed 
blazars displaying intraday variability in the y-ray band, reach- 
ing values of less than 2-3 hours (e.g. Tavecchio et al. 2010, 
Foschini et al. 2011a,b, Vovk & Neronov 2013). Specifically, it 
is interesting the case of the flat-spectrum radio quasar (FSRQ) 
PKS 1222-1-216 (z - 0.432) with the minimum observed time 



scale of 1 .0+0.2 hours (~ 0.7 hours intrinsic) measured by Fermi 
during the outburst of 2010 April 30 (Foschini et al. 2011b). 
After about 1.5 months, on 2010 June 17, the MAGIC Cerenkov 
telescope detected the same source at hundreds of GeV with a 
variability time scale of 10 minutes (Aleksic et al. 2011). The 
two outbursts were also marked by a spectral change: during the 
April one, the y-ray spectrum was soft, with a cut-off and the 
energy of the highest detected photon equal to ~ 23 GeV; in 
the second outburst, on June, the spectrum was hard and extend- 
ing up to a few hundreds of GeV. This has been interpreted as a 
change of the location of the dissipation: within the broad-line 
region (BLR) during the April outburst and outside it during the 
June one (Foschini et al. 201 lb; but see Tavecchio et al. 201 1 for 
another inteipretation). 

Recently, Ghisellini et al. (2013) reported about the case of 
PMN J2345-1555 (z = 0.621): for the first time, the multiwave- 
length coverage allowed to understand that the hardening of the 
y-ray spectrum was linked to a shift of the synchrotron peak at 
higher energies in the X-rays, confirming the change of its nature 
from "red" to "blue" (see Ghisellini & Tavecchio 2008 for an ex- 
planation of blazar colors). Other hints of changes in the blazar 
color were reported for PKS 2155-304 (z = 0.116) in 2006 
November (Foschini et al. 2008, but no y-ray data were avail- 
able) and PKS 1510-089 (z = 0.36) in 2009 March (Abdo et al. 
2010, D'Ammando et al. 2011). In the latter case, the peak of 
the synchrotron shifted from infrared to optical/UV frequencies 
and the source was detected at hundreds of GeV by HESSfQ), 
although the X-ray spectrum did not change (Abdo et al. 2010). 

Yet PKS 1510-089 displayed a very interesting activity on 
201 1 October-November with y-ray flux measured by Fermi ex- 
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Fig.l. Global light curve of PKS 1510-089 in the 0.1 < £ < 100 GeV energy band with 1-day time bin. Time starts on 2009 
January 1 (MJD 54832.0). Vertical dotted lines indicated the 1 January of each year Vertical dashed lines delimited the period of 
the three outbursts analyzed in the present work. Outburst A: 2009 March; Outburst B: 20 11 October/November; Outburst C: 2012 
February /March. 



ceeding 10"^ ph cm"^ s ' and the shorter time scale of variability 
of the order of one hour (Brown 2013, Saito et al. 2013). There 
was no multiwavelength coverage, because the source was ap- 
parently too close to the Sun, except for radio observations re- 
porting an extraordinary increase of the flux density at different 
frequencies (Nestoras et al. 201 1, Orienti et al. 201 1, Beaklini et 
al. 2011). Nevertheless, it is possible to draw useful inferences 
from the study of Fermi LAT data alone. Here we propose a re- 
analysis of the y-ray emission in the same period, but with the 
addition of a different method to extract the light curves with 
time bins smaller than three hours. This allowed us to investigate 
time variability on scales shorter than those reported by Brown 
(2013) and Saito et al. (2013). A comparison with other out- 
bursts (2009 March, 2012 February /March) suggests also some 
interesting features of the jet of PKS 15 10-089. 



PKS 1510-089 has redshift z = 0.36 (Thompson et al. 
1990). In the usual ACDM cosmology with the latest measured 
values of the Hubble-Lemaitre constant Hq - 70 km s"' Mpc"' 
and Oa = 0.73 (Komatsu et al. 201 1), the luminosity distance is 
then 1934 Mpc and 1 arcsec corresponds to 5.1 kpc. 



2. Data Analysis and Discussion 

Fig. [1] shows the y-ray emission as measured by ferm//LAT 
of PKS 1510-089 over a period of about four years. In the 
present work, we focused mostly on the outburst B (2011 
October/November), but we compared it also with the outbursts 
A and C (2009 March and 2012 February /March, respectively) 
in order to draw useful inferences on the source properties. 

The analysis of the Fermi/LAI data has been performed 
following the standard procedure^ for light curves with time 
bins greater than 3 hours (i.e. two orbits), which is the min- 
imum integration time to have a complete coverage of all the 
sky (Fig.[T]i. In addition, we also built light curves with smaller 
time bins (Fig. |2]|4| using the method described in Foschini 
et al. (2011a,b). In that case, the time bin is equal to one 
Good-Time Interval (GTI), which in turn is less than one or- 
bit (~ 95 min) and can have variable duration depending on 
the pointing direction. In both cases, the flux is estimated by 
using the gtlike task of the LAT Science Tools. The meth- 
ods were applied by using the most recent version of the soft- 
ware and calibration: LAT Science Tools v. 9.27.1, in- 
strument response function (IRF) P7SOURCE_V6, isotropic back- 
ground isop7v6source.txt, and Galactic diffuse background 



E.g. |http://fenm.gsfc.nasa.gov/ssc/data/analysis/scitools/| 
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gal2yearp7v6v0 .fits. PKS 15 10-089 and the contaminating 
sources within 10° from its sky position were modelled by using 
a simple power law F{E) oc E^^, where F is the photon index. 
The threshold of Test Statistic (TS, see Mattox et al. 1997) for a 
meaningful detection was set to 9 (~ 3cr). 

The peak flux was measured on 201 1 October 19 19:06 UTC 
(MJD= 55853.79595, source ontime ~ 2.2 ks) as (4.4 + 0.3) x 
10-5 ph cm-2 s-i, with T = 1.98 + 0.04 (TS = 598) and 119 
counts recorded. This corresponds to a luminosity of about 2 x 
10'*'' erg s^' making PKS 1510-089 one of the brightest blazar 
ever detected at y rays (the record remains to 3C 454.3, with 
3x10^" erg s"' on 2010 November 20; see Foschini et al. 2011a). 

The light curves were scanned searching for the minimum 
time of doubling/halving flux between two consecutive points: 



F(to) 



- 2-('i-'o)/'" 



(1) 



where F(?o) and F(ti) are the fluxes at the two consecutive times 
to and f 1 , respectively, and t is the extrapolated doubling/havling 
time scale. To have a meaningful measurement, we have set a 3cr 
threshold in the flux change between the two consecutive mea- 
surements. The results are displayed in Table [1] and are limited 
to tiie shortest values of the order of one hour or less. Table [1] 
reports also the values of t as measured from the light curve 
with 3-hour time bin (not shown), to compare and confirm the 
results pubHshed by Brown (2013) and Saito et al. (2013). We 
noted that the shortest r are for increasing flux, while the small- 
est value for the decay time is an upper limit of < 5 hours. This 
is consistent with the results of Nalewajko (2013), who found 
that PKS 1510-089 has the most skewed flares, with short rise 
time and long decays. 

The shortest doubling time is ~ 20 minutes (Tint = 20 + 
1 min), measured during the outburst B (Fig. |3) between 
MJD= 55852.52144 (2011 October 18 12:30 UTC) and MJD= 
55852.59130 (2011 October 18 14:11 UTC). This corresponds 
to a size of the emitting region equal to r ~ 4 x 10'"* cm (by 
assuming 6 = 10). This value has to be compared with the gravi- 
tational radius of PKS 1510-089, which is r, ~ 3 x 10'^^ cm (by 
assuming a mass of the central black hole of 2 x 10** Mq, from 
Xie et al. 2005). 

During that episode, PKS 1510-089 jumped by almost two 
orders of magnitude in less than one orbit, from a flux of (0.75 + 
0.66)x lO-*" ph cm-2 s"' (TS = 22) to (31 +9)x 10"'' ph cm^^ s"' 
(TS =95). Interestingly, the photon index at low flux was hard 
(F = 1.2 + 0.4) and became softer at high flux (F = 1.9 + 0.2), 
although the highest detected photon energy was ~ 14 GeV. No 
photon above 100 GeV was detected during the entire studied 
period of the outburst B to confirm a really hard spectrum ex- 
tending to hundreds of GeV or with a cut-off at tens of GeV, as 
expected if the dissipation occurs inside the broad-line region. 
Neither Cerenkov observation nor X-ray follow-up (to search for 
a synchrotron peak shift) was possible at that epoch, because of 
the Sun constraints. 

To better assess if the lack of very high energy photons could 
be due to the very short exposures and small LAT effective area 
above 100 GeV, we compared with the 2009 March outburst 
(outburst A, fig. |2]i. This event was already described in other 
works (e.g. Abdo et al. 2010, Marscher et al. 2010, Tavecchio et 
al. 2010). The authors reported always a rather soft photon index 
of the y-ray spectrum (F > 2). Abdo et al. (2010) reports a hard- 
ening trend of the photon index as the flux increases, but still 
from ~ 2.6 to ~ 2.2. Therefore, comparing these information 
with the present analysis (specifically. Fig. |3] bottom panel), it 
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Fig. 2. Outburst A (2009 March): Light curve in the 0.1 < E < 
100 GeV energy band with GTl-bin. (top panel) Flux in units of 
[10"^ ph cm"^ s"'] (bottom panel) Photon index. 
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Fig. 3. Outburst B (2011 October/November): Light curves of 
PKS 1510-089 in the 0.1 < £ < 100 GeV energy band with 
GTl-bin. (top panel) Flux in units of [10^^ ph cm"^ s"'] (bottom 
panel) Photon index. 
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Fig. 4. Outburst C (2012 February-March): Light curve in the 
0.1 < E < 100 GeV energy band with GTl-bin. (top panel) Flux 
in units of [10^^ ph cm"^ s '] (bottom panel) Photon index. 



seems that the outburst A (2009) was different from the outburst 
B (201 1). However, in 2009 March (A), the blazar was detected 
at hundreds of GeV by HESS ('), thus requiring a hard pho- 
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Table 1. Summary of doubling/halving times measured. Date is in YYYY/MM/DD; to and t[ are in [MJD], while half-bin times are 
in hours; fluxes are in units of [10^^ ph cm"^ s"'] and the TS of the detection is indicated between parentheses; the significance S 
of the flux difference is in [cr]; the absolute values of the observed characteristic time scale r (positive for decay, negative for rise) 
and the intrinsic one Tint = t/(1 + z) are in hours. 



date 


to 


bin 


ti 


bin F{to) F(ti) 


S 


T 


"^int 


Outburst A (2009) - GTI time bin 


2009/03/31 


54921.18046 


0.79 


54921.24677 


0.79 2.44 ± 0.53 (34) 7.88 + 1.04(18) 


5.3 


-1.36 ±0.71 


-1.00 ±0.52 


Outburst B (2011) - 3-hr time bin 


2011/10/18 
2011/10/19 


55852.43752 
55853.68752 


1.5 
1.5 


55852.56252 
55853.81252 


1.5 0.66 ± 0.58 (22) 11.19 ± 3.10 (60) 
1 .5 8.52 ± 5.95 (37) 44.30 ± 4.73 (609) 


3.4 
6.0 


-1.06 ±0.09 
-1.82 ±0.21 


-0.78 ± 0.07 
-1.34 ±0.15 


Outburst B (2011) - GTI time bin 


2011/10/18 
2011/11/01 
2011/11/04 


55852.52144 
55866.73018 
55869.71600 


0.65 
0.67 
0.77 


55852.59130 
55866.79917 
55869.78141 


0.61 0.75 ± 0.66 (22) 30.93 ± 9.33 (95) 
0.79 1.06 ±0.75 (38) 5.67 + 1.18(35) 
0.79 1.85 ±1.25 (15) 8.12 ± 2.06 (51) 


3.2 
3.9 
3.0 


-0.45 ± 0.03 
-0.99 ± 0.02 
-1.06 ±0.05 


-0.33 ± 0.02 
-0.73 ± 0.01 
-0.78 ± 0.04 


Outburst C (2012) - GTI time bin 


2012/03/03 


55989.90896 


0.63 


55989.97872 


0.65 0.76 ±0.47 (12) 4.88 ± 0.49 (72) 


8.5 


-0.90 ±0.11 


-0.66 ± 0.08 



ton index at MeV-GeV energies. Therefore, we have reanalysed 
the LAT data starting from 2009 March 20 (MJD 54910.0) and 
ending after a couple of weeks by using the GTI-bin light curve 
(Fig. |2]i. We noted that the shorter doubling time is ~ 1 hour 
and the photon index became now sometimes harder (T < 2), as 
expected. 

The blazar was detected at hundreds of GeV also by MAGIC 
during the outburst C (2012 February/March, Lindfors et al. 
2013). Therefore, we have reanalysed also the LAT data for 
the period from 2012 January 29 00:00 UTC (MJD 55955.0) to 
April 3 00:00 UTC (MJD 56020.0). The light curve is displayed 
in Fig. |4] The shorter time scale is < 1 hour (Table [U and was 
measured while changing from a low hard (F = 1 .5+0.3) to high 
soft (F = 2.2 + 0.1) flux, just as happened in 2011. It is worth 
noting that this event occurred during one successful exposure 
with MAGIC (cf Fig. 5 of Lindfors et al. 2013). 

Therefore, given the similarities of the three outbursts, we 
concluded that: (i) it is reasonable to think that PKS 1510-089 
could have been detected at hundreds of GeV also in 201 1 (i.e. 
the photon index as measured from the LAT data analysis are 
really hard and the lack of very high energy photons is due to 
the short exposures); (ii) the y-ray spectrum can be hard on very 
short time periods (< 3 hours), while when integrated on longer 
exposures, the soft component dominates. 

3. Final Remarks 

The present giant outburst of PKS 15 10-089 seems to be differ- 
ent from that of PKS 1222+216 (cf. Foschini et al. 201 lb). The 
latter displayed in 2010 two giant outbursts: the first one with a 
cutoff at tens of GeV and doubling flux time scale of ~ 1 hour; 
the second one extending without changes in the spectral slope 
up to hundreds of GeV and variability of ~ 10 minutes as de- 
tected by MAGIC (Aleksic et al. 201 1). The two outbursts were 
separated by about 1.5 months, more than sufficient to think to a 
drift of the dissipation region from inside to outside the BLR. 

Instead, the 2011 outburst (B) of PKS 1510-089 reported 
in the present work (see also Brown 2013, Saito et al. 2013) 
displayed different properties. The 0.1 - 100 GeV flux well ex- 
ceeded the value of 10"^ ph cm"^ s"' three times during the pe- 
riod of one month, with very short doubling times recorded in the 
first two bursts and hard photon indexes (Fig. |3)- Particularly, 
the first two episodes, which were separated by about two 
weeks, displayed similar spectral and variability characteristics. 
Therefore, it is unlikely that this was due to a drift of the dis- 



sipation region. Either a stable multizone structure, where the 
outburst are due to disturbances moving down the jet (Marscher 
et al. 2010), or a structured jet, with compact blobs responsible 
of the shortest and most powerful flares (Tavecchio et al. 201 1), 
are a more plausible explanation. 

Despite the limitations of the present work due to the lack 
of MW coverage, one firm conclusion can be, and has to be, 
stressed. This episode, compared to other similar cases, shows 
that there could be significant differences from source to source, 
and even from event to event of the same source. Therefore, in- 
dividual case studies are still a powerful drain of information to 
understand the nature of relativistic jets. It is of great importance 
to report as many as possible of these events, in order to set up 
a significant archive. In this scenario, the continuous monitoring 
of the y-ray sky performed by Fermi is a priceless font of data. 
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